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Part II: Optimization of electrolysis cells, using an electrolyzer 
for peroxodisulfuric acid as an example. 

Communication from the VEB Eilenburger Chemie Werk, concern of the 
Combine VEB Chemische Werke Buna, management area Research and 
Development, and of the Martin-Luther University Halle-Wittenberg, 
departments for Mathematics and Chemistry. 

The calculation of electrochemical reactors on the basis of balance 
equation systems [1] is still limited at present to simple cell 
geometries. In addition it requires the knowledge of a large num- 
ber of material- or system-related parameters which, in real cases, 
are generally unknown and which can only be determined with great 
experimental expenditures. On the other hand, in the application 
of combined statistically and physically based mathematical models, 
results can be obtained in a short time with relatively small ex- 
perimental expenditures for the determination of the empirical 
constants introduced into the model, which can make an important 
contribution toward optimization and perfection of electrolyzers 
of complicated engineering construction. 

Such models were developed in the VEB Eilenburgwer Chemie Werk and 
have been used for some years successfully for the optimization 
of engineering and economic characteristic values [2], We now plan 
to present here the model of an electrolyzer for peroxodisulfuric 
acid and to analyze it mathematically. Its application for engi- 
neering and economic optimization will be investigated in detail. 


* Numbers refer to foreign pagination. 
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Process-theoretical model 


In the modelling of the electrolyzers we shall dispense with a de- 
scription of the cell design as well as the processes in the interior 
of the cell; for this we refer to part I [ 3 ]. We use only relations 
between the external process variables to set up the model/ whereby 
the effect of thermal-engineering parameters is based on physical 
causes. The variables# which can be taken into account to set up the 
models for a fixed cell design and constant electrolye composition# 
are the current strength I# the dosing speed and the external 
temperatures# i.e.# the cooling water temperature and the ambient 
temperature #2 . In the subject application case for the production 

of hydrogen peroxide it makes sense to vary the dosing rate and the 
current strength always in the same proportion to each other# so 
that the final concentration of the peroxodisulfuric acid falls 
into the range most favorable for the subsequent hydrolysis stage. 
Thus the dosing rate is not an independent variable; the effect of 
the current strength then contains the effect of the changing flow 
velocity and the average dwell time. The relations between the re- 
maining independent variables and the quantities# which describe 
the electrolysis process# such as electrolysis temperature 0 ,cell 
voltage U# current yield n and specific electric energy con- 
sumption E g are then given sufficiently accurately by the 
simplified interrelation model shown in figure 1 . 

To set up the mathematical model it is now necessary to describe 
each interrelation by a function# using data obtained by measure- 
ments. In order to keep interferences with the operational electro- 
lysis process as small as possible# we only made measurements here 
which for fixed external temperatures and fl 2 (0) for the different 

values of current strength I determined the associated values 
for the electrolyte temperature ,r , the cell voltage U and the 
current yield . In place of the current strength I we shall 
use as variable in the following the current density i referen- 
ced to a reference area F (e.g.# the cathode area). One can now 
represent the desired target quantities in a purely statistical 
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way as functions of the starting variables. It turns out, however, 
that it is exactly the stepwise buildup of the model suggested by 
the interrelations with the use of some physically justified 
relations, which leads to an essentially more informative model. 

Initially the current yield >1 and the cell voltage U are intro- 
duced as functions of the current density i for fixed external 
temperatures ^ (0> and 

>1 =- i- + .v 2 «» i -!■ « 8 <0), U = a 0 i + b 0 ( 1 ) 

whereby the coefficients »•„«», a», fro are determined from the measured 
values by the method of least squares. The cell temperature is 
determined with neglect of reversible thermal changes from a 
simplified energy balance 

y.i( ft - d{) + x 2 (ft ~ ft 2 ) = Fi[U - - E 2 °) - i 

- E 2 °] ( 2 ) 

in which and *2 denote the heat transfer coefficients, referenced 
to one electrolysis cell, in the cooling fluid or in the environment 
and Ei° = 2,2 v, AY* -- 1,2 v the standard EMK of the anodic peroxide or 
oxygen formation. To make possible a simple solution for , 0 , we con 

sider >; here as approximately constant with a value of 0.7 and use 
for U the relation from (1), whereby the temperature dependence- 
of U is thus neglected. Under these assumptions it is possible 
to determine the values *i (0) and *2 (0) for the electrolysis cell con- 
sidered with the method of least squares from the values measured 
for ft . 


The solution of ( 2 ) for 0 produces 


F i (ao i 4- 60 — 1 >9) 


* 1 + *2 


- + ft 3 nut ft 3 = 


XI ftl + X2 fts 
xi + xs 


( 3 ) 


whereby the designation is used as abbreviation for the average 

external temperature. In order to now obtain the current yield and 
the cell voltage U as functions of the two variables i and , 
the expressions 


V =Mi) +/ 2 ( 0 ), U = ? 1 (t) + g-,(ft) (4) 


are set up with separated variables. The summands depending on 
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are determined by integration from the approximately known 
temperature coefficients 


8 n /89 = —0,0006 9, 8U/89 = —0,01 (5) 

The summands depending on i can be determined in that one 
demands the agreement of (4) and (1)/ if one inserts in (4) the 
external temperatures !V 0) and the heat transfer coefficients y i i0> 
and used as basis in formulae (1). If in the thus obtained 

functions for rj and 0 the temperature & is substituted as 
per (3), one finally obtains n and O as functions of the 
starting variables in the form 

*1 — *3 + x 2 * + 4" a '4 i 3 + 

U = ci 2 -\-ai-+-b (6) 

with the coefficients 

Xl = a;iO» + a (1 - D 2 ) + 2 a A (9 3 <°> - D d 3 ), 
x 2 = a - 2 («) + 2 a B (0 3 «» - D 9 3 ), 

x 3 = a' 3 <°> + a (0 3 <°> 3 - # 3 2 ), (7) 

= 2 AB — Z> 2 ), a 5 = a A 1 (1 — IP), <% = 0,0003 


and 

c = 0*4(1— D), a = a 0 + 0 B(l - D), 
b = bo + 0 ( 0 3 < o) - 03), 0 =-= 0,01 ( 8 ) 

with the abbreviations 

E «o 0’(6o - 1,9) xi«» + *:>«» 

- «i<°> + x/o) ’ h - x X m 4 - x 2 1 °> ’ D ~ 

If and U are known# one obtains the specific electric— energy 
consumption as 

E, = c 0 u/t) mit co = 2 1\,/M = 270,26 A • li/kg (9) 


For better visualization of the formulae obtained we have shown in 
figure 2 the current yield and the specific electric-energy con- 
sumption as a function of the current density for different values 
of the cooling-water temperature 9i for 9 2 20 "o . In addition we 

have plotted the current yield- current density functions for con- 
stant electrolyte temperatures 0 . The plotted curves are based on 

values which are encountered in a commercial electrolyzer for 
peroxodisul furic acid and namely 
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0!<«> 
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= 3 V, 

f/ 2 (0) r _, *20 (\ .a<*» 
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(10) 


Engineering-economic optimization model 

The optimum operating point of a facility is determined not only by 
the engineering-physical data, but also to a considerable extent by 
the economic data. In order to set up the associated mathematical 
model, an analysis of all the expense types affecting all the total 
costs of an electrolysis process is necessary; these must be pre- 
pared in the form of cost functions. 

Among the fixed costs, which are independent of the current density 
and the amounts produced we count the basic equipment expenditures 
for the electrolyzers and secondary facilities in the form of annual 
amortizations and the fixed operating costs, such as for example 
wages and salaries. The variable costs include; 

- the costs for the electrolysis current; they are proportional to 
the power consumption and are calculated with an electric-power 
cost d^ and the efficiency d 2 of the DC current 

generation as 

K\ — ci </i Ui/<h with = S <>00 h/u 

- the costs for general repairs whose frequency increases pro- 
portionally to the electrolysis current 

K 2 — </;{ i 

- the energy costs for the cooling facility which are proportional 
to the quantity of heat which is picked up by the cooling medium 
along a square meter reference area F 


K 3 -■ Cl 1/1 '1.1 y.i (/> •— Ui)/r 

(here d^ is the ratio of the expended electric power to the 
quantity of heat picked up and removed by the cooling fluid). 


5 


/681 


ORIGINAL PAAA 13 
OF POOR QUALITY 


- the costs proportional to the annual production quantity m 


/' i — ■ 


With rs ” c i/ c * ~ 3', 13 kj» «iA 


Thus the total costs encountered per year per square meter of 
electrolyzer area are 


K A() -*■ <-*1 *t\ i f*| '/] '/[ (/7 - ■ j ) ^ /' j- 

! ca«/3< »/ (li) 


The analysis of the total costs for a realistic engineering 
electrolyzer showed that the costs for the electric 
power with almost 60% represent the greatest share, followed by 
the production-proportional costs with almost 25%. The target 
functions important for an optimization of the electrolysis 
process now are the specific costs k per unit quantity pro- 
duced 

'• ' K k/o-j. > v) n-) t ! 


and the operating result per year and per square meter of elec- 
trode area 

HE <1$ m — K /i 


whereby dg denotes the calculated price for the electro- 
lysis product. 


For better visualization all target functions are shown in figure 3 
as functions of the current density. Here those values were used 
as bases for the economic parameters, which can be encountered in 
a production persulfuric acid electrolyzer, and namely 

A'o = 1 000 M/a- 111 2 , dx — 0,1 5 M/k\V • h, 

< I 3 0,3 M/a- A 

dx = 0,9, d t = 0,4, <?5 = 0,1 M/kg, </« = 0,6 M/kg 

One can see in the figure that the various target functions always 
have an optimal value which, however, is assumed for different 
current density values. 
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Evaluation of the model 


It was shown that that the derived model reflects the electrolysis 
process in the range, interesting from a production standpoint, 
between 900 and 1500 A/m with sufficient accuracy. It should 
not be extrapolated beyond this range since for lower current den- 
sities, e.g., the linear expression (1) for the cell voltage is no 
longer correct and since for higher current densities, e.g., the 
expression (6) for the current yield at constant electrolyte tem- 
perature decreases again while from a physical-chemical standpoint 
one could expect that it approaches a limiting value. 

Between the optimum current density values for the here considered 
target functions of the process-theoretical and production-economic 
optimization, there exist a few general relations which are always 
valid independently of the individual parameter values. They either 
satisfy the inequality chain 

( ? 'opt ) b , < ( ! opt)ij £» (*op t)t < (iopt)]}E < (*opt)m (14) 

or the inequalities 

(*opt)£ # < (*opt)r; < (»opl)m, (topt)fc tS (('upt)i;, 

(*"opl)fc < (*o <C (£op»)»i ( 15 ) 


We shall dispense here with proofs for all inequalities; we shall 

only examine in detail the in practice particularly interesting 

inequality (toptk- < (*opt)«£. From the necessary extreme-value requirement for 

* * / 
k = K:m there follows K m - K m = 0. If one solves this for m 

and inserts this into the derivative of BE, one obtains 

Ii'm K’ 

BE' = (/(j»r — K' = d o —j£ K' = — (tin m — K) — 

= -- BE > 0 (1(5) 

If we assume that the operating result BE is positive, which 
only seems to make sense here, then BE is positive at (*opt)* 
i.e., the operating result is a function still growing there. 
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Thus if in an arbitrary electrolysis facility one is operating at 
a minimum cost setting/ then the operating result can always be still 
increased. Reversely the setting for maximum yield is not the most 
cost-effective one. Therefore a new facility should be laid out such 
that the projected production quantity is reached at a cost minimum/ 
but that a variation of the current density is possible in the range 
between the extreme values for both target quantities. 

The optimum current density values/ whose general setting with 
respect to one another has just been characterized/ can be easily 
calculated by an approximation method in every real case via the 
necessary extreme value conditions. Figure 4 shows for the parameter 
values/ used as example here/ the optimum current densities and the 
relative function values at the optimum, referenced to # 3 = 8 "C , 
as a function of tf 3 ,or d x < 

After the optimum for a given parameter value has become known, it 
is often interesting to find out how the optimum changes when the 
parameter value is changed. Thus one must calculate the change of 
the implicitly determined optimum current density as the result of 
the necessary extreme value requirement and of the extreme value of 
the function. Since only small changes of the parameter are considered, 
so as not to go beyond ther validity range of the developed model, 
the development of a power series provides a useful approximation 
solution for this purpose. Thus we calculated in this way for the 
concrete case of a few parameter values the relative changes of the 
optimum current density listed in table 1 or the optimum target /g 

function value for a 1% change of the corresponding value of the 
parameter. One can see that a change of the starting parameters 
affects the operational result considerably more than the specific 
costs. Furthermore it becomes clear that the location of the opti- 
mum current density values depends greatly on the economic charac- 
teristic values, such as the price of electric power and the price 
of the product. It can also be seen that a change of the statis- 
tically determined coefficients xi <0) , x 3 <®> of the current yield- 

current density function substantially affects the location and 
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the magnitude of the optimal values. Thus the dependence of these 
coefficients on design quantities is to be investigated in the 
continuation of the present project. 

Now special interest must still be given to the dependence of the 
target functions on the physical-engineering parameters found in 
the model. The effect of the cooling-water temperature on some 

target functions/ as can be seen from figure 4, shows that a low 
cooling-water temperature is the most favorable even if the 
cooling costs are taken into account. Since limits are set for a 
lowering of the cooling-water temperature by the performance of 
the cooling facility/ other measures for lowering the electrolyte 
temperature are important. Therefore we shall next consider the 
efect of the heat transfer coefficient *1 . An improvement of the 

heat transfer properties increases the costs for the construction 
of the electrolysis cell; as an example let us assume that the 
fixed costs can be determined by the formula 

A'o = A'oo + rfo (y-i - y-i m ) . 

For do = 1/ 2/ 3 M . K/a . m^ . W we calculated with this 
assumption the minimum specific costs and the cost-optimal current 
density as shown in figure 5. One can see that the minimum costs 
first decrease with increasing heat transfer coefficient xi , but 
then increase again. The thus determined best value for the heat 
transfer coefficient provides information for an optimum design 
electrolysis cell from a thermal-engineering standpoint. 

One can obtain an improvement of the heat transfer properties 
especially easily by an enlargement of the cathode surface which 
essentially coincides with that of the cooling area. If one intro- 
duces an enlargement factor p into the model and if one takes 
into account the dependence of the fixed costs of p in an 
expression corresponding to formula (17), then the calculation 
of the specific costs as a function of p provides results which 
are similar to those shown in figure 5. From this one obtains in- 
formation about an optimally designed electrolysis cell with 
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respect to the cathode area. 


The results and conclusions derived from the model show that even 
with a relatively simple model and with relatively low measurement 
expenditures one can derive important information not only for 
optimization of the operating conditions of an existing electrolyzer, 
but also for continued design developments and that here an analysis 
of the formally mathematical structure can already provide infor- 
mation about which changes of the target functions are at all 
compatible with the design of the model. The here presented model 
makes possible with continuing scientific understanding of the 
electrolysis process a stepwise continued development in that 
further interrelation expressions are explained physically and that 
the effect of engineering design parameters is revealed by the 
model. In this way the extrapolation capability can be improved 
and an optimization of the electrolyzer with respect to additional 
parameters becomes possible. 


Summary 

A combined statistically and physically justified model is pre- 
sented and evaluated with an example of an electrolyzer for 
peroxodisulfuric acid. An engineeering economic optimization is 
made with the aid of an analysis of all costs affecting the 
electrolysis process. By means of a mathematical analysis of the 
model we drew generally applicable conclusions concerning the 
change in position of the optimum with respect to the various 
investigated target functions as the result of changes of the 
individual design-caused and economic parameters. 
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symbols 

a o h o 

a,b,c 

BE 

cj 

dj 

E s 

Ej° 

F 

F o 

I 

j 

K 

K j 

k 

M 

m 

U 

xj (o) 

r j 

0 

l>i w) 

Di 

n 

xj (o) 
xj 


explanation 

coefficients of the cell voltage function (starting cell) 
coefficients of the cell voltage function (arbitrary cell) 
operating result per cell and per year 
conversion factors (constants) 
cost factors 

specific electric power consumption 
standard EMK 

reference area of one cell 
Faraday constant 
current strength 

current density (referenced to the cell area F) 

total costs per cell and year 

cost shares per cell and year 

specific costs 

molar mass 

production quantity per cell and year 
cell voltage 

coefficients of the current yield functions of the starting cell 

coefficients of the current yield functions of an arbitrary cell 

electrolyte temperature 

external temperatures (starting cell) 

external temperatures (arbitrary cell) 

current yield 

heat transfer coefficients (starting cell) 
heat transfer coefficients (arbitrary cell) 
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In part I of this series, which was published in Chem. Techn. 29(1977) 148 
to 154, the part of the table belonging to figure 1 was accidentally not printed 
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Figure 1 Interrelation model 

1- cooling water temperature ambient temperature $ 

2 - current strength I (current density i) 

3- bath temperature 4-current yield 5- bath voltage 
6- spec, electric power consumption E s 



Figure 2 Relation between current yield and specific electric power 
consumption and the current density for different cooling water 
temperatures 0i <»*«>) .. »o '*c) . 

1- current yield functions 

2- spec, electric power consumption 
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Figure 3 Current density dependence of the investigated target 
functions - 20*0 



Figure 4 Dependence of the relative extreme function values, 
referenced to the starting temperature »:»(«> and of 
the optimal current densities on the external temp- 
peratures and the cooling-water temperature ih 

(for 2 ~ 20*0) 
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Figure 5 Dependence of the minimal specific costs and the cost- 

optimal current densities on the heat transfer coefficient 
for cells with design changes having different fixed cost 
factors nth ■'> <■■ »v> io'u 



| *i«» 

X 


*1 

do 


K 0 

((opt)#£ 
(iopt)jfc 
U^>max 
— **mia 

— 1,059 
— O.Sll 
1 — 2,973 

: — 0,531 

1,219 

0,817 

5,135 

1,092 

0,095 

— 0,006 

0,708 

0,202 

— 0,246 
— 0,076 
— 1 ,684 
— 0,639 

0,051 

— 0,592 

— 0,237 

0,306 

3,552 

0,075 
— 0,199 
— 0,093 


Table 1 Relative change of the optimal current density and the 
optimal target function values during changes of the 
parameter value by 1%. 
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The model is presented of an electrolyzer for per- 
oxodisulfuric acid, and it is analyzed mathematically. 
Its application for engineering and economic 
optimization is investigated in detail. The math- 
ematical analysis leads to conclusions concerning 
the change in position of the optimum, with respect to 
the various target functions due to changes of the 
individual design-caused and economic parameters. 
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